COMMENT ⊗   VALID 00002 PAGES
C REC  PAGE   DESCRIPTION
C00001 00001
C00002 00002		SUBROUTINE ELLIPSE(A,B,PSI)
C00003 ENDMK
C⊗;
	SUBROUTINE ELLIPSE(A,B,PSI)
	N=IFIX(A*50.)+12
	IF(B.GT.A)N=IFIX(B*50.)+12
	PHI=(6.2831853/360.)*PSI
	X=A
	Y=0.0
	THETA=6.2831853/FLOAT(N)
	C1=COS(THETA)
	C4=A/B
	C6=SIN(THETA)
	C2=C6*C4
	C3=C6/C4
	C=COS(PHI)
	S=SIN(PHI)
1	CALL PLOT(A*C,A*S,3)
	DO 10  I=1,N
	U=X*C1-Y*C2
	V=X*C3+Y*C1
	X=U
	Y=V
	X1=X*C-Y*S
	Y1=X*S+Y*C
2	CALL PLOT(X1,Y1,2)
 10	CONTINUE
	RETURN
	END